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Abstract: 

^ Coulomb effects on tt and K spectra in relativistic nuclear collisions are investigated. At collision energies around 
■ 1 GeV the ratio of 7r~ to Tr~^ is enhanced several times at low transverse momenta but less at ultrarelativistic en- 
' ergies. We describe the ratios at SIS, AGS and SPS energies with simple analytic models as well as more elaborate 
7-H numerical models incorporating the expansion dynamics. The Coulomb effect depends on the properties of the source 
Q.^ , after the violent collision phase and provides information on source sizes, freeze-out times, and expansion velocities. 
(— I Comparison with results from HBT analyses are made. Predictions for tt^ and at RHIC and LHC energies are 
given. 

> 

' PACS numbers: 25.75. +r, 25.70.Pq 



I. INTRODUCTION 

An asymmetry in the number of opposite charge pions has been observed at intermediate energies Q and has 
recently also been identified in heavy ion collisions at energies around 1 A-GeV H,^, 11.4 y4-GeV and 158 A-GeV 
1^. The effect is strongest at lower beam energies. The ratio of negative to positive pions at low pion momenta is 
~ 3 at SIS energies, but only ~ 1.6 at AGS and SPS energies for the central collisions of heavy nuclei as Au or Pb. 
A possible cause for this effect is the Coulomb interaction between the produced pions and the positive charge from 
the reaction partners. 

Earlier work |^-^ has described the tt^ /n^ enhancement/reduction by Coulomb attraction/repulsion from a 
static source. In this paper we underline the importance of the source dynamics and of the longitudinal and transverse 
expansion of the fireball. We emphasize the different reaction kinematics and dynamics that evolve from 1 A-GeV 
to 160 A-GeV. At the lower collision energies nuclear matter is fully stopped and expands relatively slowly in all 



directions whereas at higher energies much of the initial motion survives and the system expands faster preferentially 
in longitudinal direction. 

Since the Coulomb forces influence the matter essentially after freeze-out, the asymmetry in the number of 
differently charged particles can be directly related to the freeze-out parameters of the matter. Therefore, a detailed 
description of the violent initial part of the collision is unnecessary, and it is sufficient to consider the state of the 
nuclear matter at the time Tf of freeze-out. Similar to particle interferometry measurements the observation of 
Coulomb effects can help to determine the dynamics of heavy ion collisions. 

We will develop approximate models by assuming simple geometrical shapes of the charge distribution at freeze- 
out. In these models we can estimate analytically the ratio of tt^ to 7r+ as a function of the transverse momentum. 
In an earlier publication ||lo[] we have presented a model of an expanding pion gas created in heavy ion collisions 
at ultrarelativistic energies. This model is well applicable in a wide range of collision energies. We will use it for 
calculating particle ratios as a check of the simpler models. 

The manuscript is organized as follows. In section II we study slowly expanding and spherically symmetric sources 
from which pions escape rapidly. Such a scenario may be applicable to heavy ion collisions around 1 A-GeV and we 
compare to SIS data. In section III, we study rapidly expanding systems from nuclear collisions at ultrarelativistic 
energies, which may apply to collisions at the AGS and SPS. We compare to recent AGS and SPS data. In section IV 
we investigate the situation for particles at very low momentum. In section V we briefly overview the main ingredients 
of the dynamical model developed in ref. |lO) and compare. In section VI we compare freeze-out radii to those obtained 
from HBT analyses. In section VII we discuss the influence of particles created at late stages from resonance decays, 
and in section VIII we consider / ratios. Finally, predictions for charged pion and kaon ratios at RHIC and 
LHC energies are given. 



II. SLOWLY EXPANDING CHARGE 



In nuclear collisions charged particles feel a Coulomb force due to the net positive charge originating form the 
protons in the colliding nuclei. In particular the light produced particles as tt^, A"^, ... moving slowly with respect 
to the net positive charge are significantly affected by the Coulomb field. We will concentrate on pions, which are 
affected strongest by the Coulomb field since they are the lightest particles, but all results apply to any particle of 
mass m as, e.g., kaons or protons. 

The net charge of the source resides mainly in protons also after the collisions. Expansion of the "cloud of charge" 
can thus be considered slow when typical proton velocities Vp are much smaller than pion velocities w^r . At SIS energies 
the energies of protons from central collisions and at midrapidity follow approximately a Boltzmann distribution with 
temperature Tp ~ 100 MeV Only slow pions with v^^ < Vp are affected while for pions with kinetic energiesfl 




*We will in the following use the convention Ti = c = = 1- 



the expansion can be ignored. At SIS energies this hmit is thus ~ 20 MeV and Eq. (|l]) apphes to aU SIS data shown 
in Fig. |l|. Thus the charge expansion can be ignored for the SIS data. 

Let us therefore consider a particle {t:^,K^, ...) with mass m and energy Eq = ^ w? + which is formed in 
a spherically symmetric source with net charge Z evenly distributed over the volume of radius K. We will for the 
moment assume that the charge is static and will later extend the analysis to an expanding charge. The Coulomb 
potential varies from Ze^ /R at the surface to 3/2 times that value at the center with the average 

The energy and momenta of tt'^'s are unaffected by the Coulomb field and define the reference values, Eq and poj 
whereas the tt^ are accelerated/decelerated in the field. As the volume expands slowly with respect to the fast pions, 
the final tt^ energy is 

E = Eq±V,. (3) 

The effect on the particle distribution functions can be derived from particle conservation, i.e., the number of 
particles before Coulomb effects are not changed after, dNo=dN. The tt^ pion distributions at freeze-out d^No/d^po 
now have to be transformed by changing their energy and momentum following (^). The final distribution is 

m ^ d^fpo ^ (mPoEo ^^-^ 
d'^p d^po d^p d^pQ p E 

The original dN'^ /d^po can be taken from tt" data around 1 GeV. These pion spectra are to a good approximation 

given by dN^/d'^po ^ exp{~Eo/T) with slope parameter T. With final tt^ energy and momentum given by and 

using for both tt^ and 7r+, one obtains the ratio (see also ||8|,|l^) 



TT- dN-ld^p ,TT- , ( 2K\ S + K \{E + YcY -m^ 

-' Uxp -— ——J (5) 



TT+~dN+/d^p V+' 'V T J E-Vc\] {E-Vc, 
The square root denotes the ratio of the initial ir" and 7r+ momenta when they are detected with final energy E. 
It leads to diverging rate as the kinetic energy {E — m) approaches Vc, 7r~/7r+ cx {E-rn- VcY^I'^. The average 
(7r^/7r+) is the ratio of the total number of negative and positive pions produced. At collisions energies of 1 GeV 
most pions are produced through the A resonance, the ratio is 

— = — ^ ~ 1.94 , 6 

where iV and Z are the neutron and proton numbers of the Au + Au nuclei before collision. 

The simple formula (^) gives a good description of the tt^/tt^ ratio for the 1 GeV Au + Au collisions at SIS 
and allows us to determine both T ~ 75 MeV and Vc — 27 ± 2 MeV, see Fig. |[ These two parameters can also 
be estimated independently. The measured single particle slope T ~ 75 MeV agrees nicely with our estimate. 
According to Eq. the value of Vc — 27MeV and Z = 2 x 79 x 0.75 (14 % centrality) corresponds to an average 
source size of i? ~ 8.0 fm, which is larger than the geometrical radius of the initial participant zone of the two colliding 
Au nuclei. 

The corresponding baryon density at pion freeze-out is ^ 2A/ (AttR^ /3) ~ 0.85/Oo, i.e., somewhat smaller than 
normal nuclear matter density po- We will see that this result agrees with the investigations in the dynamical model 
described in section V. In the analysis of this data within the transport model [|l3| a similar value of the Coulomb 
energy was obtained, although at smaller density of ~ 0.6po- 



III. RAPIDLY LONGITUDINALLY EXPANDING SYSTEMS 



At the higher AGS, SPS, RHIC and LHC energies the systems expands rapidly in particular in the longitudi- 
nal direction along the beam axis. For systems expanding more rapidly than the pion velocity the net charge is 
reduced with increasing time resulting in decreased Coulomb effects. At ultrarelativistic speeds one also needs to take 
retardation effects of the Coulomb fields into account. 

Our starting point is the retarded time component of the electromagnetic potential from a charge Q moving with 
velocity v along the beam axis jlj] 

'^{r^,t)^A'^ ^ (7) 

y'v^t-' + (1 — v')rj_ 

where r± is the distance from the beam axis and the time is chosen such that t = when the charge is closest to 
the point r±_ (see Fig. The vector part of the potential points along v and is therefore longitudinal and does not 
contribute to the transverse electric field that affects transverse momenta. The longitudinal component even vanishes 
for an approximately constant rapidity distribution since the contributions of negative and positive velocities cancel 
each other. 

In ultrarelativistic heavy ion collisions the charged particles have large longitudinal but small transverse momenta. 
The distribution of longitudinal velocities v can be obtained from the measured rapidity distributions since v = tanh(y) 

dv 1 — dy 

In central Ph + Ph collisions at SPS energies the net proton rapidity distribution is rather constant at midrapidities, 
dN'^P^dy = 37±2 

Assuming an approximately flat charge rapidity distribution, the electric field becomes 

dN''' r^dv 
d^(4 + i;2(<2_r2^))3/2 

(9) 



dN"'' 2 r I 



dy t 

In relativistic heavy ion collisions the upper and lower velocities range from ±c .0 In the Bjorken scaling picture IT^ ] 
t counts time from the moment t = of the total overlap of the collision partners. 

As can be seen from (|^) the electric potential is dominated by the slow charges (v ^ 0) simply because they 
provide a Coulomb field for a longer time. Thus pions with rapidity y mainly feel the Coulomb field from charged 
particles at similar rapidities. If one instead assumes a charge particle rapidity distribution peaked at target and 
projectile rapidities, (i.e. the initial charge distributions of the colliding nuclei) then one finds that the electric 
potential is suppressed by a factor 7c„is ''^ 100 near midrapidities at SPS energies as compared to the constant 



^The exact expression is obtained by replacing the time t in (j^ by v^t'^ + r\{l ~ v^), where vo is the lower or upper velocity. 
Here, vo < c limits the longitudinal extent of the charge cylinder and thus removes the t — singularity for nonvanishing rx . 
Since the times and distances relevant for Coulomb effects are t c± R/c and ^ R respectively, where R is the size of the 
nuclei, we can safely approximate iiq — c = 1. 




rapidity distribution. In contrast, assuming a static central charge the push given to the positively charged pions is 
greatly overestimated. 

As above we assume that the pions are trapped in the hot matter until freeze-out at time Tf . After this time the 
pions can propagate freely only influenced by the Coulomb field. During the collision the particle rapidity distributions 
change due to interactions. However, at late times the retarded potential (^) arises mainly from the free streaming 
phase or the flowing phase in the period between the collision and freeze-out. According to Gauss law the electric field 
in the interior is reduced by an area factor r\/E?p where Rf now is the transverse size of the system at freeze-out, 
i.e., 

^ATch o I r\ < Rf I 

(10) 

A charged pion of transverse momentum p_l,o at freeze-out and final momentum p± receives a momentum change 
or "Coulomb kick" 

/•oo 

Ap± = p± - p±,o = ±e / E{r^,t)dt, (11) 

where the ± refers to positive and negative pions, respectively. If the Coulomb kick is small so that the pion velocity 
changes little, the position of a tt* at freeze-out ry with velocity developes according to 

r^Vf + v^it-Tf). (12) 

For an ensemble of pions we obtain the momentum change by averaging over the freeze-out position rf in Eq. ([fl]). 
The result will generally depend on both Rf and TfVj_. We find, however, that the momentum change is only a slowly 
decreasing function of the freeze-out time. For a freeze-out time of r/ — 8fm/c, Rf ^ lOfm and with pion velocities 
as given by the low m± data, we obtain the estimate 

p,= Api ~2e2— — — , 13 
dy Rf 

in the direction of ±px. The variation with freeze-out time is ±10% for vj_Tf /Rf in the range 0.5 — 0.8. 

It is important to note that the difference between results for the static charge of Eqs. (||J^) and the rapidly 
expanding charge of Eqs. { pUjl^ ). In the former case the field decreases at large distances as whereas in the 
latter it decreases as t^^rj^^. As a consequence, the energy is shifted by an amount proportional to the charge Z in 
the former case while in the latter the momentum is shifted by an amount proportional to dN'^^ /dy. In both cases, 
however, the scale of the shifts is determined by the source size Rf at freeze-out. 

As above we can derive the Coulomb effect on the transverse particle distribution functions 

dN _ dNp pog 
d'^p± d2po,_L P± 

The pion spectra are well reproduced by dN/d'^p± oc exp{—m±/T) with T ~ 150 MeV at both AGS and SPS energies 



Defining mj^ — y/ m? + {p± ± pc)'^ and using Eqs. (11) and ( |l3| 

_dN^/d'^p± /""^x f — p± + p, 



7r+ dN+ /d'^p± 7r+ \ T J P± — Pc 



exp^— ^ i:^-^. (15) 



The pion ratio predicted by (|T5|) is compared to AGS and SPS data in Figs. || and ^. Here the experimental values 
(7r^/7r+) ~ 1.2 and 1.05 at AGS ||l^ and SPS energies, respectively, have been used. These values are closer to 
unity than at SIS energies due to the abundant pion production not connected to A-resonances. The experimental 
data is well reproduced with Ap^'~^^' ~ 20 MeV/c and Ap^^^ ~ 12 MeV/c. From the measured proton rapidity 



distribution, dN^^g/ dy ~ 70 and dNgp^/dy ~ 37, and Eq.(13) we can now extract the size of the systems at freeze- 
out: R^'^^ ~ 10 fm and Rf^^ — 9 fm. The sizes are substantially larger than the geometrical radii Rgeom — 6 ^ Rpb 
fm of the collision zones (for 14% centrality) indicating that significant expansion takes place before freeze-out. 

For the 5+5* collision at the SPS dN^ /dy ~ 4. Assuming i?^ ~ 4 fm we find a rather small value for pc — 4MeV/c, 
and the corresponding Coulomb effect is thus very small as seen in Fig. 

IV. SLOW PARTICLES 

Slow pions, i.e. pions with velocities Vtt smaller than the expansion velocity of the net charge which is typically 
given by the proton velocities Vp, are less affected by the Coulomb field. According to Eq. (|l|) such slow pions should 
have kinetic energies below 20 MeV at SIS energies and transverse kinetic energies below 40 MeV at SPS energies. 

We can estimate the effect of expansion by simply assuming that the protons solely determine the expansion of the 
charge. If we assume that the proton velocities are distributed according to a nonrelativistic Boltzmann distribution 
with temperature Tp, the charge should therefore be reduced by a factor 

^e.p ^ Cexp(-lmX/rp) v^-'dv il exp(-x2), n - 2 1 ^^^^ 

exp{-^mpV^/Tp) v"-^dv | er/(x) - (2/0F)x exp(-a;2), n = 3\ 



where x = \J ^mpv'^/Tp and erf(x) is the error function. The exponent n = 2 is for cylindrical geometry with 
transverse expansion and n = 3 is for spherical expansion. Reducing the charge (or proton rapidity density) by this 
factor reduces the Coulomb kick Vc (or pc) is diminished by an amount which now depends on the particle energy. 
This fact needs to be properly included in the Jacobian dpo^^/dp^ of the resulting charge particle distribution (Eq. 
(|l^)) The result of transverse expansion is shown in Fig. || for 7r^/7r+ with pc — 15 MeV/c and in Fig. ^ for / . 

The factor ( p^ ) overcorrects for transverse expansion since it completely removes the Coulomb corrections for 
Utt — 0. However, even for = the pions are still accelerated by the Coulomb field which decreases as i^^. We 
therefore consider slow pions embedded in a charge expanding with velocity Vp. Pions with velocities smaller than 
Vp stay inside the expanding source in the potential decreasing with time. As the net positive charge Z or 2dNP /dy 
expands spherically or transversely with velocity Vp, the radius increases with time 

R^Rf + Vp{t^Tf). (17) 

Here, Rf is the radius at freeze-out radius at freeze-out time tj. The pion moves according to Eq. ( |l2| ) and we will 
find that the Coulomb effect leads only to minor changes in pion momenta, hence to leading order we can assume 
that is constant. 

In the case where pions move slowly in the electric field 

E(t) = Ze^ . (18) 



of a sphere with a constant charge density the net momentum change is 



poo 

Ap = ±e / E{t)dt = ±- 

Jre 



Irs "^RfVp 



Vr, Rl 



(19) 



for positive and negative pions, respectively. Retardation effects are ignored for simphcity which is a rather good 
approximation for spherical systems where the charge moves primarily in radial direction (see also appendix A). The 
average of the freeze-out position vanishes, (rj) = 0, if there is no transverse flow. Therefore, averaging over will 
to first approximation only lead to a smearing of the 7r+ and tt^ distributions and by the same amount. When taking 
ratios this effect cancel out, and we therefore ignore it in the following. The final momentum is thus 

P = P0(1±C3), ^3 = ^^, (20) 

i.e., the positive and negative pion momenta are scaled up and down respectively. 

The case of rapid longitudinal expansion was studied in the previous section. The resulting electric field is given 
by Eq. ( p^ ) and from Eq. ( p^ ) we find the change in transverse momentum 

Ap± = p± - p±.o = ±e / E{r^,t)dt~±e''——— — . (21) 
Jrf ay Rf Vp 

The latter approximate expression holds under the condition VpTf^Rf. For small p± we can approximate pj_ — mv± 
and thus find 

p^=p^,o(l±C2) ; C2^e^— -. (22) 

ay mVpH / 

The Coulomb effect on particle distribution functions can be obtained from particle conservation, i.e. dN — dNo 
is the same for initial and final pions, 

dN dN" 

(1±C„)-". (23) 



d"-p d"po 
The ratio of negative to positive pions is then 



exp , (24) 



7r+ dN+fd'^p 'n+'yi-CnJ V T 

where = \J m? + p^(l ± C„)^. Notice that C„ <C 1 is implicitly required since we approximated the trajectory 
r(t) by assuming constant velocity in the electric field. 

Introducing the kinetic energy E-^in = E — m (for n — 3) and Ekin — fn± — m (for n — 2), the ratio decreases 
linearly as 

^-(^)(l + 2"^^«(l-^)) ■ (25) 

for small kinetic energies, E^in <C m. The ratio is thus enhanced by 2nC at small m± and decreases linearly on a 
scale determined by T. 

It is interesting to compare our approximation with the exact analytic result found by Ayala and Kapusta pl| ] 
treating the motion of nonrelativistically positively charged particles in a Coulomb field from a homogeneously charged 
expanding sphere. With their parameters {Z—158, Rf — 8 fm, Vp—0.2) we obtain C3 — 0.7 for kaons. This leads to 



a reduction of low momentum /f+ by a factor (1 + Ca)^"^ = 0.2 similar to their value of 0.15 determined from Fig. 
3a in pl| . However one should be aware that the suppression factors behave asymptotically differently in the limit 
Vp 0. 

Next we estimate the crucial parameters C2 and C3. In central Au + Au collisions at SIS energies the proton 
velocity can be estimated from the slope of proton spectra ~ STp/nip ~ 0.3. With a central charge oi Z = 122 (15% 
centrality) we obtain C3 ~ 0.3 that according to (24) leads to a large enhancement by a factor ^ 6 at zero momentum. 
Unfortunately, at SIS energies no measurements are available at very low pion energy. For SPS we assume cylindrical 
geometry n = 2. Consequently, going from C3 to C2 we replace Z by 2dNP/dy. With ~ 2Tp/mp ~ 0.6 we obtain 
C2 =0.12 and find a 7r~/7r+ ratio ( |2^ ) of 1.8 in agreement to the data. 

As mentioned the derivation applies not only to pions but to all particles, for example, kaons. As the kaon is 
heavier we expect from ( |2^ ) that the Coulomb enhancement of the /K^ ratio is smaller by a factor of mT^/mK ^ 
1/3. 

Experiments at the AGS and SPS |]l8| indicate that transverse flow is important in collisions with very heavy 
nuclei. Transverse flow decreases the m± slopes observed in m± spectra of pions, kaons and protons. The measured 
"effective temperatures" are thus larger than the intrinsic temperatures. In our analysis transverse flow has several 
effects. The temperature employed in the Boltzmann spectrum should not be the effective slope parameter, T ~ 200 
MeV but the intrinsic temperature, T ~ 140 MeV. Furthermore, flow leads to larger Vp and nonvanishing (r/) in (p^). 
Particles with transverse momentum p_L from a thermal source of temperature T with transverse flow urj_/Rf arrive 
on average from a point displaced from the center by (r/) = p±uR/2T. Transverse flow thus adds a term to ([l9| ) for 
charged pions and usually reduces C„. 

We can summarize the discussion of Coulomb effects by addressing the important scales entering. In the spher- 
ically symmetric case (n=:3) the Coulomb kick is given by Vc — ztZe^/Rf. For slow expansion the particle kinetic 
energies are shifted (see Eq. (|^)) by this amount, i.e., the relevant dimensionless quantity is Vc/mT^v'^ instead. For rapid 
expansion Vp 3> Wjr the kinetic energies are instead scaled (see Eq. ( |20| ) by an amount proportional to Vc/riiT^Vp. In the 
dimensionless quantity is replaced by Vp simply because it is the fastest particle velocity that matters. A similar 
behavior is found for the cylindrical expansion (n=2) except that in this case the relevant quantity is the momentum 
instead of kinetic energy due to the rapid longitudinal expansion and the Coulomb kick is Vc — {e'^ / Rf)dNP /dy. 
For slow expansion the momentum is shifted by Vc (see Eq. (13) and the corresponding dimensionless parameter 
is Vc/TOTrWjr- For an expansion faster than the pion velocity the momenta are scaled by an amount proportional to 
Vc/rriTrVp (see Eq. (|2^)). Again the pion velocity is replaced by the proton velocity since the fastest particles determine 
the magnitude of the Coulomb effect. 



V. DYNAMICAL MODEL 



After the previous discussion of the typical properties of the Coulomb effect in different energy regions we apply 
a more detailed model jioj which allows finer details of the freeze-out scenario. Inspired by the Bjorken picture of 
approximate boost invariance we assume a cylindrical geometry and include both longitudinal and transverse flow. 
We introduce the following parameters: the width Aych of the rapidity distribution of the charge (usually a Gaussian 



distribution is taken), the thermal equihbrium temperature T, the mean transverse flow velocity (/3), the source 
radius Rf, and the freeze-out time Tf. The initial charge distribution at freeze-out is assumed to follow a Boltzmann 
distribution which includes the different chemical potentials for positively and negatively charged particles. Further 
details are explained in ref. [p^ . 

After freeze-out the different species of particles, mainly nuclcons, pions, and kaons, are assumed to move freely. 
Their momentum distribution is given by a thermal distribution superimposed on a hydrodynamic flow. The influence 
of the electromagnetic forces on the emitted particles will increase the initial charge asymmetry. This influence is 
treated as a first order perturbation, assuming that the electromagnetic field is generated by the charged particles 
moving with their unperturbed initial velocities. This method is justified as long as the electromagnetic potential is 
small compared to the mean kinetic energy of the particles. 

To calculate the electromagnetic potential A'' at the space-time ~ (t, r) 

the knowledge of the currents j^{t,r) during the whole collision history is needed. Retardation effects are essential, 
especially at ultra-relativistic energies where the longitudinal motion is not fully decelerated. To describe the situation 
before freeze-out we have in ref. |l0[| used a simple hydrodynamical scenario to describe the situation before freeze-out. 
We assume a linearly increasing longitudinal and transverse flow of the matter between the overlap time t — and 
the freeze-out time t = Tf, i.e., the radius at freeze-out is Rf — Rgeom + 0.5Tf(3s, where Rgeom is the initial or 
geometrical radius of the overlap zone at the collision time and Ps the velocity at the surface of the source. The linear 
motion of the particles allows one to calculate the net Coulomb mean field as a superposition of the electromagnetic 
fields of the individual charged particles. If the particles denoted by i with charge g^, mass m.^, and four- velocity 



— (1, v)/\/l — start at space-time yf we get 

Since the calculations of retardation are very cumbersome we also discuss approximations in appendix A. 
The motion of the particles is treated relativistically in the potential (|27 



m.-^u, =q,u,^—~—y (28) 

where tj denotes the proper time of the particle. The tt^ and 7r+ spectra are finally calculated by a sampling the final 
momenta over a set of solutions obtained with different initial conditions obtained from the freeze-out parametrization. 

In Figs. we show the calculations for the three above mentioned experiments. We concentrate on the 
sensitivity to variations of the freeze-out time. The flow velocity and temperature were determined from flts to the 
m± slope of available spectra for particles with different masses, see refs. pl] , p^ , p^ . For a systematic overview we 
refer to ref. ||. 

Fig. shows the result for SIS energies which has to be compared to the schematic model of Fig. H. For a 
temperature of 75 MeV and mean transverse velocity of {(3) — 0.3 corresponding to a surface velocity of (3s =0.5 good 
agreement to the data is obtained. This gives a freeze-out radius of Rf = 6/m + Q.hTfPs— 8 fm leading to a break-up 
density of 0.8 po- This value is just the same as that obtained in section II, but larger than the value obtained in the 



transport model p3[ |. The final production time in that model is rather long ~ 25 fm/c which is partially due to long 
lived A resonances. 

Figs. ^ and|^ address the ultra- relativistic regime. Here, the measured widths of the particle rapidity distributions, 
Aj/c/i, is about 0.7 and 1.3 for AGS and SPS energies respectively. A mean flow velocity of (/3) =0.42 or (3s — 0.62 was 
used together with a moderate temperature of 140 MeV and 120 MeV, respectively. The results at AGS energies are 
relatively insensitive to the freeze-out time since the charge diffuses slower in the longitudinal direction than for SPS 
energies. A freeze-out time of 10 fm/c seems to be compatible with the data, while for SPS energies a time of about 7 
fm/c is required. We also show the relatively modest enhancement for the reaction of S* on at SPS energies, where 
the rapidity density of the charge is slightly depleted at midrapidity. For a more detailed discussion at SPS energies 
we refer to ref. |p^ . 



VI. COMPARISON TO HBT ANALYSES 

The source for pions, kaons and other particles created in high energy nuclear collisions at freeze-out is also 
studied in two-particle correlation functions. Such Hanbury-Brown and Twiss (HBT) analysis provides complementary 
information on freeze-out time and source size. Pion correlation functions are, however, affected by resonance decays 
which actually produce the majority of the pions. Resonances lead to an increase of the measured HBT radii, 
effective source sizes and life-times 24 . A significant fraction / of the pions stem from long lived resonances 
in ultrarelativistic collisions which reduces the correlation function by a factor X^"^ ^ (1 — f)"^. Experimentally, 
^TTTT ^ Q g _ Q g gj^^ so / ~ 25%. Also, long lived resonances reduce the Coulomb effect by a minor amount 
as discussed below. Though these and other effects as Coulomb corrections, detector cuts, etc., are important for a 
detailed analysis, we will ignore them for our qualitative comparison. 

The transverse HBT radii are usually referred to as the outward and sideward radii in the direction along and 
perpendicular to the transverse momentum of the two correlated pions, respectively. The outward HBT radius contains 
a number of fluctuations as duration of emission and resonances life-times and are furthermore affected by opacities 
and moving emission surfaces |p6|| . The sideward HBT radius Rg is more directly connected to the geometric size of 
the source at freeze-out. For a source of constant pion density the sharp cut-of radius is given by i?/ — 2Rs, and with 
~ 5 fm ||2^ in central Ph + Pb collisions at 160 A-GeV we obtain Rf ~ 10 fm. This number is larger than the 
geometrical radius of a Ph nucleus indicating that expansion takes place before freeze-out. 

In a longitudinally expanding system with Bjorken scaling the freeze-out time is related to the longitudinal HBT 



radius by t/ ~ Ri ^Jm^/T. From the NA44 data Ri ~ 5.5 fm |g5[ we estimate Tf ~ 8 fm/c. For comparison RQMD 
calculations predict substantially larger values for the freeze-out time, Tf — 15 fm/c [p7| . Assuming that the source 
starts out with the mean geometrical radius Rgeom — 6fm (for the 15% centrality) at the collision and develops 
transverse flow j3s up to freeze-out, we find a radius of Rf = Rgeom + Ps^f /2 — 9 fm, which is slightly smaller than 
the freeze-out radius obtain from the sideward HBT radius. We conclude that both freeze-out times extracted from 
the longitudinal and transverse HBT radii are compatible at SPS energies with the one we find from tt^/tt^ data. 
The Coulomb freeze-out radius is also comparable to the HBT radii found at AGS energies as seen from Fig. |l^. At 
SIS energies we extract the HBT radius from the average of the 7r+7r+ and tt^tt^ HBT radii by multiplying by 



the geometrical factor, Rf — \/5/3Rrms — 15fm. This radius is considerably larger than the one extracted from the 
Coulomb effects. 



VII. RESONANCE DECAY 



In ultrarclativistic collisions resonance decay into pions is important. If pions stem from resonances with a life 
time larger than the freeze-out time, their Coulomb push is reduced. This is important for the long lived resonances 
uj,ri,ri' , Kg, etc.. As mentioned above they amount to roughly / = 25% if the reduction of the HBT correlation 
function is solely caused by long lived resonances. If a fraction / is unaffected by the Coulomb field only the fraction 
(1 — /) of the distribution function in Eq. ( |l^ should be transformed. The resulting 7r^/7r+ ratio is therefore 

/T 

(29) 



The effect of long lived resonances are shown in Fig. || for / = 0.25 and pc ~ 15 MeV/c. As they reduce the Coulomb 
effect, the best fit is now obtained for a slighly larger value of Pc — 15 MeV/c. If the long lived resonances produce 
pions at different pj^ one should use a different effective m_L slope for those in (^9|). 

Resonances may decay differently into tt^ and 7r+ depending on pion momenta and thus provide another mecha- 
nism for changing pion distributions at low p^. This is the case for A and A particles. In central Ph -I- Ph collision at 
SPS energies the number of A's produced per unit rapidity at central rapidies is dN^^^ /dy = 19 ± 7 ||l^. Of these 
64% decay into tt^. Compared to dN^ /dy = 160 ± 10 the A's thus contribute ^ 8% to the tt" yield. Since the decay 
energy of pions amounts to only 35 MeV in the A rest frame they would lead to an enhancement of the 7r~/7r+ ratio 
of ~ 1.4 if included. The ratio would fall off slowly as a function of the transverse momentum. However, only at most 
a few percent of those pions can enter the detector in the NA44 measurements because the A's travel ^ 60cm before 
decaying due to the long life time. Therefore they do not contribute to the measured tt^ yield. Both the Pb + Ph 
data and the S + S data are compatible with the Coulomb effect but cannot be explained by resonances alone. The 
same conclusion was drawn from RQMD analyses . 



VIII. THE K- /K+ RATIO 



The K~ jK^ ratio provides another way to test the Coulomb effect and obtain information on freeze-out. We 
estimate Apj^ by taking the best fit to the 7r^/7r+ data, however, correcting it for the slightly smaller HBT radii for 
kaons than pions. We estimate Ap^ ^ 15 MeV/c. Furthermore we take the kaon temperature of Tk = 230 MeV from 
Ipsi and {K^ / K^) = 0.68 from We can then predict the jK^ as shown in Fig. |[ As the kaons are heavier 
than the pions their kinetic energies are less affected by the momentum change. Also, the kaons move slower and 
are more affected by transverse expansion of the net (proton) charge (see appendix A). A smaller number of protons 
remains to give the kaons a Coulomb kick and the / ratio is closer to unity as seen in Fig. ^. 

HBT analyses show that the kaon correlation function is not reduced by long lived resonances, since ~ 1. 
On the other hand, the absorption of on nucleons by K^N An is much larger than for K^. Therefore the 
initial ratio / may differ in very heavy ion collisions from light ones since the stopping is larger. The measured 



spectrum ||3^ in fact shows a decrease at low p±^ in heavy ion coUisions as compared to hght ones by almost 
a factor of two. Thus, the resulting j ratio could be even an increasing function of the transverse energy. 
Recent NA44 measurements do not show a significant dependence on the transverse energy. If the enhancement of 
the / ratio predicted in Fig. ^ should exceed the measured one, this would indicate new physics. 

Kaons may also experience a strong interaction with the mean field from hadrons. Early reports on low pj^ 
enhancement for at AGS energies - referred to as "cold kaons" - led to suggestions of a strong attractive potential 
that dominated the repulsive Coulomb potential. The cold kaons have, however, not been confirmed in later AGS 
experiments. In nuclear matter the interaction is repulsive jSSt but one cannot rule out that it becomes attractive 
in hot and dense matter. 



IX. SUMMARY 

Coulomb effects are very sensitive to the expansion dynamics of high energy nuclear collisions. At lower energies 
the expansion of the charge cloud, which consists mainly of protons, is slow compared to typical pion velocities. The 
pions experience the full Coulomb kick and the 7r^/7r+ ratio becomes large at small pion momenta as described in 
Eq. (§). At higher collision energies the expansion of nucleons is faster and exceeds that of the very slow pions. The 
Coulomb field, which the pions feel, vanishes rapidly within the typical expansion time. The resulting -k" jix^ ratio 
is only enhanced by ~ 60% at the smallest measured transverse momentum at AGS and SPS energies and can be 
approximately described by Eq. ( p^ ) . At ultrarelativistic energies the rapid longitudinal expansion and retardation of 
electromagnetic fields are important and the charge rapidity density replaces the central charge. The Coulomb field 
affects small pion transverse momenta on a scale given by the total Coulomb energy. The Coulomb effect on kaons is 
smaller than for pions due to their larger mass. 

The Coulomb effect provides information on the sources produced in relativistic heavy ion collisions at freeze-out 
independent of HBT analyses and flow studies in particle spectra. HBT methods analyse the source at freeze-out, and 
Coulomb effects occur after freeze-out. The measured enhancement of the 7r^/7r+ ratio provides a severe constraint 
on the values found by various methods. From the investigations presented here the source sizes and freeze-out 
times obtained from HBT analyses seem to be compatible with those obtained from Coulomb effects. Calculations 
using transport models like BUU or RQMD tend to find larger freeze-out times indicating the need for a consistent 
description of the whole process to calculate charge asymmetry and HBT correlations. 

At RHIC and LHC energies smaller stopping is expected than at SPS and AGS energies, and the charge density 
dN'P / dy will be smaller at midrapidity but larger near target and projectile rapidities. Given the amount of stopping 
we can try to predict the resulting 7r^/7r+ ratios at RHIC and LHC. The total {tt^ /n^) production should be very 
close to unity due to the abundant pion production. The transverse size of the system at freeze-out can be estimated 
from the HBT radii assuming that the radii continue to increase slowly with collision energy as is found experimentally 
up to energies I60A-GeV. By extrapolating we estimate i? ~ 15 fm. If we, for example, assume dNP/dy = 25 wc find 
Pc — 5 MeV/c. The resulting 7r^/7r+ ratio is shown in Fig. ^ (not corrected for resonances and transverse proton 
expansion). However, if a phase transition to a quark-giuon plasma occurs, these predictions may change drastically. 
The freeze-out may occur much later so that the source is more longitudinally and transversely expanded leading to 



a reduced Coulomb effect. 

We acknowledge very informative discussions with A. Hansen and Nu Xu from the NA44 collaboration. Support 
from The Danish Natural Science Research Council is appreciated. 



APPENDIX A: Effects of Retardation 

In ref. approximations were applied by using the Coulomb potential cj) as the only part of the four-potential. 
This potential 4> is calculated as being A'^ of (|26|), but neglecting the retardation. Thus, the equation of motion 
simplifies to 

"^^^^'^'^^N^^' fc= 1,2,3. (30) 



We demonstrate the deviations one obtains when applying the approximative method (^0|) where retardation and 
the spatial components are neglected. 

At moderate energies where the thermal velocities of the nucleons are not large this approximations works well. 
In Fig. |l| we show the effect on the tt /7r+ ratio in the transverse direction kinetic energy for a reaction of gold on 
gold nuclei. As an example we use the following parameters: Aj/ch = 0.4, T — 75 MeV, {[3) — 0.3, and r/ = 7 fm/c 
which give in Fig. |^ reasonable agreement with the data. One recognizes that the approximate treatment (^0|) is well 
applicable and leads roughly to a 10% underestimation of the ratio at an lab-angle of 44° which corresponds nearly 
to the transverse direction in the cm. system. 

Larger effects are expected at ultrarelativistic energies. In Fig. |l2| we show the calculated charge ratio at small 
momentum as a function of rapidity for a lead on lead collisions at 158 A-GeV. The solid line is calculated using 
Eq. ( p^ ) and shows a /n^ ratio at zero transverse momentum which follows nicely the rapidity distribution of the 
net charge as it is expected from our considerations in section III. The approximate method shows quite a different 
behavior with a maximum off-midrapidity. This is because the potential (f) is not able to cope with the translational 
invariance which is hidden in the correct potential (^6|). Since the calculation is carried out in the cm-system it fails 
completely for larger rapidities. In the correct treatment the action of the component A'^ is compensated by the 
longitudinal part of the magnetic field. The dotted curve shows the result if one treats the Coulomb potential as a 
genuine scalar potential which means dropping the term on the right hand side of Eq. ( |30| ) . This method is quite 
an insufficient approximation which even is not applicable for small energies as can be seen in Fig. |ll|. 
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FIG. 1. Experimental tt^/tt^ ratios for the reactions Au + Au at 1.0 AGeV as a function of the pion kinetic energy 
around 44 degrees from forward direction in the lab system (corresponding to ca. 90 degrees in cms) . Curves show calculations 
according to Eq. (jH]) for the case of a slow expanding spherical source using various Coulomb energies Vc- Using a slope 
parameter of 75 MeV a value of Vc ~ 27 MeV gives good agreement with data (full squares) [0|. 




FIG. 2. Sketch of the ultrarelativistic colhsion where the pion at the transverse distance r_L is influenced by a charge Q 
moving with velocity v in the expanding matter. 




m,-m (MeV) 

FIG. 3. Experimental -k' /■k'^ ratios Q for the reactions Au + Au at 11.6 AGeV as a function of transverse mass compared 
to calculations employing Eq. (|l3) using a constant momentum shift pa- 
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FIG. 4. Experimental tt^/tt^ ratios at SPS energies for the reaction Pb + Pb as a function of transverse mass compared 
to calculations using Eq. ( |l^ ) using a constant momentum shift pc- 




m,-m (MeV) 

FIG. 5. The effect of resonances and transverse expansion on tt^/tt^ ratios Eq. (113) (see Fig. and text). 
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FIG. 6. The predicted K /K^ ratios at SPS energies for the reaction Pb + Pb as a function of transverse mass from Eq. 
with Apx = 15 MeV/c (full curve). Dashed curve includes transverse expansion of Eq. jl^). 
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FIG. 7. Calculations of -r^ /tx^ ratios with the dynamical model for the reactions Au + Au at 1.0 AGeV as a function of 
pion kinetic energy in the center-of-mass system compared to experiment Q. Good agreement is obtained with r = 10/m/c 
for a moderate mean flow velocity (/?) = 0.3 c or r = 7/m/c for (/3) = 0.4 c. 
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FIG. 8. Experimental tt /tt"*" ratios ^ for the reactions Au + Au at 11.6 AGeV as a function of transverse mass compared 
to calculations using the dynamical Coulomb model for different freeze-out times r/. 



158AGeV Pb+Pb 
m m NA44 data 

T=120 MeV, Ay 



<I3> = 0,42 



1 .3 



'/,=3 fm/c 
T,=7 fm/c 
T,= 1 5fm/c 




0.8 



1.3 - 



1.2 



200AGeV S+S 

■ ■ NA44 daia 

T=140 MeV, Ay = 1 .3 
<^> = 0.3 

T,= 3 fm/c 

T,= 7 fm/c 

T,= 1 5 fm/c 




0.9 - 



100 200 

™±"™jt (MeV) 



300 



400 



FIG. 9. Experimental tt /tt^ ratios [g| at SPS energies for the reactions Pb + Pb and S + S" as a function of transverse 
mass compared to calculations employing the dynamical Coulomb model with three different values for the freeze-out time. 
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FIG. 10. The freeze-out radii for central Au + Au or Pb + Ph collisions from HBT and Coulomb effects. Tlie HBT data 
are average values of 7r+7r+ and 7r"7r" sideward HBT radii taken from Refs. Q, (29| and |2|,|l|j3| at SIS, AGS and SPS 
energies respectively. The radii extracted from Coulomb effects are compatible with those from HBT except at SIS energies. 
The geometrical radii of the overlap zones (for ~ 14% centrality) at time of collision are smaller indicating that significant 
expansion takes place between collision and freeze-out. 
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FIG. 11. -k' /■k'^ ratios for Au on Au collisions at 1 AGeV as a function of the center-of mass kinetic energy calculated with 
the correct retarded electromagnetic potential (solid line) in comparison with an approximative treatment retaining only the 
time component without retardation (dashed line) and using a simple scalar Coulomb potential (dotted line). 
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FIG. 12. TT /7r+ ratios for Pb on Pb collisions at 158 AGeV at transverse kinetic energies smaller than 10 MeV calculated 
with the correct potential and two approximations (see caption of Fig. ^l|) . 



